Dynamics of fractal dimension during phase ordering of a geometrical multifractal 
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A simple multifractal coarsening model is suggested that can explain the observed dynamical 
behavior of the fractal dimension in a wide range of coarsening fractal systems. It is assumed that 
the minority phase (an ensemble of droplets) at t — represents a non-uniform recursive fractal set, 
and that this set is a geometrical multifractal characterized by a /(a)-curve. It is assumed that the 
droplets shrink according to their size and preserving their ordering. It is shown that at early times 
the Hausdorff dimension does not change with time, whereas at late times its dynamics follow the 
f(ct) curve. This is illustrated by a special case of a two-scale Cantor dust. The results are then 
generalized to a wider range of coarsening mechanisms. 

PACS numbers: 64.60. Ak, 61.43.Hv 
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Fractal growth phenomena have been under extensive 
investigation during the past two decades [|]-U ■ The in- 
verse process of fractal coarsening occurs in many phys- 
ical systems. It has been discussed in the context of sin- 
tering of fractal matter || . Coarsening of fractal clusters 
by surface tension in the bulk-diffusion-controlled 
interface-controlled and edge-diffusion-controlled 
systems has been investigated. Additional examples in- 
clude thermal relaxation of rough grain boundaries |{J 
and smoothing of fractal polymer structure in the process 
of polymer collapse M] . Two-dimensional fractal finger- 
ing, observed in a Hele-Shaw cell with radial symmetry 
(for a review see Ref. pj|), exhibits coarsening at a late 
stage of the experiment. All these systems are quite dif- 
ferent, as they involve non-conserved or conserved order 
parameters, different transport mechanisms, etc. 

A crucial issue related to any phase ordering process 
is the presence or absence of dynamical scale invariance 
(DSI) |T3|. DSI assumes that there is a single dynamical 
length scale X(t) such that the coarsening system looks 
(statistically) invariant in time when lengths are scaled 
by A(i). Does a fractal cluster or a fractal interface ex- 
hibit DSI (on a shrinking interval of distances) in the pro- 
cess of coarsening? Early scenarios of fractal coarsening 
in systems with non-conserved |13 and conserved B or- 
der parameter did rely upon the hypothesis of DSI. How- 
ever, numerical simulations showed that the DSI breaks 
down during the coarsening of fractal clusters in edge- ||] 
and bulk-diffusion-controlled 0] systems. On the other 
hand, recent simulations of smoothing of a fractal poly- 
mer during collapse |10|, and of interface-controlled frac- 
tal coarsening under a global conservation law fl4{ , do 
support DSI. Therefore, a question arises about possible 
universality classes of fractal coarsening. 

Even if DSI holds, the fractal dimension may or may 
not change with time. Early fractal coarsening scenar- 
ios |l3|,^| assumed that it remains constant (again, on 
a shrinking interval of distances). Experiments on sin- 
tering of silica aerogels (a convenient way of investigat- 
ing fractal coarsening) have been inconclusive. Some 



of them y gave evidence in favor of constancy of the 
fractal dimension during coarsening, while others |]l5| re- 
ported a significant change of the fractal dimension with 
time. Another evidence for a significant decrease of frac- 
tal dimension with time was found in experiments on 
thermal annealing of ferroelectric thin films of lead zir- 
conate titanate fl6|j . In this experiment, the fractal di- 
mension remained constant at early times, and decreased 
to its final value at intermediate times. Numerical sim- 
ulations of a variety of coarsening systems with differ- 
ent growth laws showed that the fractal dimension does 
not change with time. These simulations include bulk- 
diffusion-controlled |||(|, edge-diffusion-controlled ||[§, 
and interface-controlled systems. 

It is remarkable that in so many systems with widely 
different coarsening mechanisms the fractal dimension re- 
mains constant during the dynamics. Therefore, one is 
tempted to look for a general scenario that would ex- 
plain this fact and that would be insensitive to specific 
coarsening mechanisms. The simple multifractal coars- 
ening model developed in this paper has this property. 
In addition, this model is the first attempt to address 
the multifractal properties of fractal coarsening. 

We shall consider a very simple model of a coarsening 
fractal system. In this model, the initial condition for 
the minority phase is an ensemble of droplets that rep- 
resents a geometrical multifractal. We will then assume 
that the smaller droplets shrink and disappear indepen- 
dently, according to their sizes, and consider discrete time 
dynamics. Using a well-known theorem of multifractal 
geometry, we will establish the dynamical behavior of 
the Hausdorff dimension of this simple coarsening sys- 
tem. This result will be illustrated in a special case, 
when the droplets are distributed in the form of a two- 
scale Cantor dust Jj],||]. Employing the size distribution 
function of this fractal set [jlTj, we will follow the dy- 
namical behavior of the rf-measure in two characteristic 
limiting cases and show that the Hausdorff dimension's 
dynamics in this example are consistent with the general 
result. Then we will relax the discrete time assumption. 
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Furthermore, we will show that the results are essentially 
independent of the details of the coarsening dynamics as 
long as the minority-phase droplets do not merge or break 
up. 

The minority phase of our model represents, at zero 
time, a big but finite ensemble of droplets that form a 
non-uniform recursive fractal Q with a constant density 
distribution in the iS-dimensional space. Let us index the 
droplets in the m-th generation of the fractal according 
to their radii. Thus, all the droplets with index k have ra- 
dius R m (k) and form a subset of the whole fractal which 
we denote by S m (k). The smallest droplets have index 
k = and radius R m (0), which is the lower cutoff of the 
fractal. The largest droplets have index k = m and ra- 
dius R m {m), which is the upper cutoff. One can work 
with a size distribution function n m {k), which is simply 
the number of droplets with radius R m (k), and use it to 
compute the Hausdorff dimension of the fractal (see Ref. 
fl7|| , where this was done for a two-scale Cantor dust). 

Any non-uniform recursive fractal with a constant den- 
sity distribution can be described as multifractal in ge- 
ometrical sense (see Ref. [pi, p. 66). In this case one 
can introduce the measure of the subset S m (k) in the 
following way: 



Hm{k) 



^k=o nm (k)Rm{k) 



(1) 



where R m (k) are the radii of the droplets divided by the 
size of the system. The Holder exponent of the elements 
of the subset S m (k) is defined by fLSj] 



a m {k) = 



In Hm{k) 
In Rm(k) 



(2) 



The f{a) curve for the fractal is constructed in the fol- 
lowing way 



/(«) 



In n m (k) 
In R m {k) 



(1 < k < m) , 



(3) 



where k is supposed to be expressed through a with the 
help of the equation a m (k) = a. (We assume that this 
equation gives a one-to-one correspondence between a 
and k.) /(a) is assumed to have a single maximum which 
is attained for a = ao, so that /(ao) is the Hausdorff 
dimension of the whole fractal. We also assume that 
f(a(k)) is the Hausdorff dimension of the subset S m (k). 
This assumption, widely used in the physical literature, 
was rigorously proved in the case of a two-scale Cantor 
dust pQ] , and also for a class of other multifractal mea- 
sures J2l| . 

We now turn to describe the dynamics. We assume 
first that the droplets shrink and disappear indepen- 
dently, according to their radius only, and also simplify 
the governing dynamics by introducing a discrete time r 
(later we will relax these two assumptions). In the first 



time step r = the smallest droplets with radius i? m (0) 
disappear, while the sizes of the other droplets do not 
change. In the next time step r = 1 the elements with 
radius i? m (l) disappear, and so on. The set of droplets 
that survive after each step of these dynamics obviously 
remains self-similar (on a shrinking interval of distances). 
The main result of this paper is the following behavior of 
the Hausdorff dimension D as a function of the discrete 
time r. For r < fc(ao) D does not change: D(t) = Dq, 
where Dq is the Hausdorff dimension of the initial con- 
dition. For r > k(a ) D(t) = f(a(k min )) where k min {r) 
is the fc-value of the smallest droplets which have not 
yet disappeared by time r. This dynamical behavior is 
illustrated in Fig. 1. 




FIG. 1. Hausdorff dimension D of the ensemble of droplets 
versus discrete time r (solid curve). The dashed curve is the 
f(a(k)) curve at r = 0. Dq is the Hausdorff dimension at 
r = 0, while ko = fc(ao) is the value of k for which f(a(k)) 
has its maximum. 

The proof of this result is based on the following the- 
orem: the Hausdorff dimension of a union of two fractal 
sets Si and S2 with fractal dimensions Ds t and D$ 2 , re- 
spectively, is D = max(Ds 1 , Ds 2 ). (See, for example, j|, 
p. 17.) In the last time step of the dynamics, r = m, 
the coarsening object consists of the subset S m {m) alone, 
and its Hausdorff dimension is /(a (to)). In the previous 
time step t = to — 1 the object consists of two sub- 
sets : S m (m) with Hausdorff dimension /(a (to)), and 
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S m (jn — 1) with Hausdorff dimension f(a(m — 1)). It 
follows from the shape of the f(a) curve of the initial 
fractal that f(a(m — 1)) > /(a (to)). Using the theorem, 
we get D(t = to — 1) = f(a(m — 1)). 

More generally, consider time r = ko + s, where s 
is a positive integer and fco = k(oto)- At this time we 
can regard the object as a union of two fractal sub- 
sets S m (ko + s) and S m (m > k > fco + s + 1). Here, 
S m (m > k > fc + s + 1) is the union of all subsets S m (k) 
with fc = s + 1, to. It is also the whole coarsening ob- 
ject at time r = fco + s + 1. Assume by induction that 
D(t = fc + s + 1) = D(S m (m > k > fc + s + 1)) = 
/(a(fco + s + 1)). It follows from the shape of the f(a) 
curve that f(a(k Q + sj) > f(a(k + s + 1)). Hence, us- 
ing the theorem, we conclude that D(t = fco + s) = 
f{a(k{) + s)). Since fco + s is the index of the smallest 
droplets which have not yet disappeared we can write 
this result as 



D(r > fc ) = f(a{k min )) . 



(4) 



The dynamical behavior of the Hausdorff dimension 
at times t < fco can be found in a similar way. For 
t = fco — 1 the object can be considered as a union of 
two fractal subsets S m {m > k > fco) and S m (ko — 1). 
It follows from Eq. (|) that D(r = fc ) = D(S m (m > 
k > fco)) = Do- From the shape of the /(a) curve we 
get D = /(a(fco)) > /(a(fc - 1)). Therefore, D(t = 
fco — 1) = Do- More generally, for any time r = fco — s 
the coarsening object can be considered as a union of 
the two fractal subsets : S m (m > k > kg — s + 1) with 
Hausdorff dimension Do and S m (ko — s) with Hausdorff 
dimension /(a(fco — s)). From the shape of the f(a) curve 
we deduce f(a(ko — s)) < Dq. Hence, by using the above 
theorem we conclude that D(t = fco — s) = Dq. More 
generally, we can write: 



D(t < fc ) = /Mfco)) = D . 



(5) 



Let us now turn to the particular case, when the en- 
semble of droplets at zero time represents a two-scale 
Cantor dust. Recall that the initiator of this fractal is 
an D-dimensional cube of unit side length. The genera- 
tor consists of m cubes of side l\ and rii cubes of side 
?2 where 1% > l\. In each step of the fractal construc- 
tion every full cube is replaced by the properly rescaled 
generator. After the last step of the construction, which 
is the rn-th step, all the cubes are replaced by spherical 
droplets with the same size as the cubes. 

Now assume that this two-scale Cantor dust under- 
goes the simple coarsening dynamics described earlier. 
For convenience, we will compute the time-dependent 
d-measure of a two-scale Cantor dust which consists of 
cubes (the ones which were replaced by the spheres after 
the TO-th generation of the construction) . The only differ- 
ence in the computed d-measure will be in a d-dependent 
prefactor. Since this prefactor is independent of fc and to, 



it will not affect the dynamical behavior of the (i-measure 
and the Hausdorff dimension. 

The d-measure of the TO-th generation of a two-scale 
Cantor dust can be written as |17| 



M d = / n m {k)R d m {k)dk = 
Jo 

1/2 r rn 

exp [g[k)] dk . 



27rfc(m — k) J J 



where 



g(k) = — fcln ( -t ) — (to — fc) In 



TOfl2^2 



mnili 



(6) 



(7) 



and R m (k) = l™~ k l\ is the size of the cubes in the subset 
S m (k). The function exp[<?(fc)] has a (sharp) maximum 
at 



fco(rf) 



mlf + n 2 li 



(8) 



For d = Do one can show that fco (-Do) = fc(«o) = fco- At 
time t = k m i n the d-measure of the object is 



M d (r)= / n m (k)Ri(k)dk. 
Jk min (r) 



(9) 



As long as k m i n (j) -C ko(d), one can apply the saddle 
point argument used in Ref. |0| and conclude that 



M d (r « fc (d)) ~ M d (r = 0) = {mlf + n 2 l d 2 ) 



(10) 



This implies that during the early stages of the dynam- 
ics the d-measure remains, with an exponential accu- 
racy, constant. Correspondingly, the Hausdorff dimen- 
sion which is computed by solving the same equation 



mlf + n 2 li = 1 



(11) 



for d, does not change with time. This is in agreement 
with Eq. (||) obtained in the general case. 

On the other hand, when fco(d) <C t = k m i n <C to, the 
behavior of M d (r) is quite different. Since for fc > fco(d) 
g(k) is a decreasing function of fc, the main contribution 
to the integral in Eq. (^) comes from a close neighbor- 
hood of fc = fc m m(r). Therefore, in Eq. (^) we can 
expand g{k) around fc = k m i n (t) to the first order and 
get 



M d 



Tl 7n {kjnin^ RrriiJ^rniri) 
| Q {kmin) | 
h(£min,d) [y(S,min,d)] r ' 



,1/2 



where £ mi „ = k mm /m 

i-e 



mm \ { s« 



n 2 l d 



(12) 



(13) 
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and 



h 1 (€min,d) = [2-!r£ min (l - imin)] 1 ' 2 I In 



(1 — £mm)™2^2 



(14) 



The Hausdorff dimension of the subset labeled by 
is given by 



(1 ^ £min) In 'l + £min m ^2 



It follows that 



Tim {kmin) 



( m ^ 



1/2 



(15) 



(16) 



Hence, we obtain the following expression for M d in the 
limit of feo(d) <C t <C to: 



m l/2 



d-/(a(fc mi „)) 



(17) 



We see that, up to logarithmic corrections resulting from 
the factor h{£, m i n ,d), the d-measure obeys a power law 
of Rm(kmin) with a time-dependent exponent. 

Eqs. (|l^)-([li|) allow one to calculate the Hausdorff 
dimension of the ensemble of droplets in the limit of 
ko(d) <C r <C m. Taking the logarithm of both sides 



of Eq. (12) and dividing by to, we get 



InM, 



I In 

to to 



,1/2 



+ Hy{Uin,d)]. (18) 



As to 3> 1, the first term on the right hand side of Eq. 
( |l8| ) can be neglected. Therefore, the Hausdorff dimen- 
sion is determined by solving the equation 



y(tmin,,d) = 1 



(19) 



for d. The solution is just the Hausdorff dimension of 
the subset £ m i„ given by Eq. (15). Therefore, D(t) = 
f(a(k min )) for fc ^ T ^ m i m agreement with the 
general result (^). 

We now show that the assumptions of a discrete time 
and of the independent shrinking of the droplets can be 
relaxed. It is sufficient to assume only that the dynam- 
ics of each droplet are determined by its radius (and 
possibly by a time-dependent "critical radius", charac- 
terizing some mean-field interaction between droplets). 
We should also assume that the droplets do not merge 
or break up. Under these assumptions the number of 
droplets in each subset is constant (until the droplets dis- 
appear) and all the droplets belonging to the same subset 
have the same (time-dependent) radius. In addition, we 



forbid nucleation, which is a standard assumption for a 
coarsening stage Q . 

Let us denote the radii of the droplets belonging to the 
fc-th subset at time t by R m (k,t). The <i-measure of the 
fc-th subset at time t is given by : 



M d {m,k,t) = n m (k)R m (k,t) . 
This can be rewritten as : 

Rm {k, t) 



M d (m,k,t) = M d (m,k,0) 



Rm(k,0) 



(20) 



(21) 



where R m (k,0) and M d {m, k,0) are the initial values of 
the radii and d-measure. Since the initial condition is a 
geometrical multifractal, M d (m, k, 0) can be expressed in 
the following manner : 



M d (m,k,0) 



Y(-,d,{Pi}) 

TO 



(22) 



where the function Y and the parameters {Pi} charac- 
terize the initial fractal condition considered. (In our 
example of the two-scale Cantor dust the role of the func- 
tion Y was played by y, while the set of parameters {Pi} 
included n\,n2,h and I2.) Substituting (22) into (pl|), 
taking the logarithm of both sides, and dividing by to we 
get 



]nM d (m,k,t) 



In 



Y{-,d,{P t }) 

TO 



d , 

+ — In 


Rm (k, t) 


TO 


Rm(k,0)_ 



(23) 



For typical coarsening mechanisms R m (k,t) grows with 
time slower than exponentially. For example, this is 
true for non-conserved dynamics (model A) and for the 
Lifshitz-Slyozov theory of conserved dynamics (model B) 
p2| . Therefore, when R m (k,t) > R m (k,0) the second 
term on the right side of Eq. ( |23| ) is negligible at to ^ 1 . 
Similarly, it is negligible when R m (k,t) < R m (k,0) 
as long as R m (k, t) is not exponentially smaller than 
R m (k,0). Eq. (E3) becomes inconvenient in the case 
of shrinking droplets at the moment of their disappear- 
ance. Eq. (j2l]) shows, however, that the d-measure of 
such droplets vanishes. Hence, the d-measure of the 
A;-th subset does not change during the coarsening dy- 
namics until the droplets belonging to this subset dis- 
appear. Consequently, the Hausdorff dimension of this 
subset does not change until its disappearance. We have 
therefore shown that the results of our simple discrete 
time coarsening model apply to a wide range of coarsen- 
ing mechanisms. It should be noticed that for a system 
with weak multifractal properties our model predicts that 
the fractal dimension remains approximately constant at 
all times. Therefore, this model provides a simple ex- 
planation to the observation that the fractal dimension 
does not change in a wide range of coarsening processes 
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In summary, we have considered a simple model of 
coarsening disconnected droplets forming a geometrical 
multifractal. We have shown that at early times the 
Hausdorff dimension of the system does not change, 
whereas at late times its dynamics follow the /(a) curve 
of the initial multifractal distribution. These results are 
insensitive to the particular coarsening mechanism. We 
hope that they will motivate experimental investigation 
of multifractal aspects of fractal coarsening. 
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Science Foundation, administered by the Israel Academy 
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